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I. INTRODUCTION 


This  report  describes  the  work  done  on  Contract  No.  N001 23-76-C-1371 
for  the  year  I August  1977  to  31  July  1978. 

The  overall  program  is  divided  into  three  areas.  These  are  1)  basic 
scattering  code  development,  2)  reflector  antenna  code  development  and 
3)  basic  theoretical  studies  to  support  the  first  two  areas.  The  following 
sections  describe  the  research  accomplished  in  each  of  the  three  areas 
mentioned  above. 

II.  PROGRAM  SCOPE 

The  scope  of  the  work  under  Contract  No.  N001 P3-76-C-1371  is  to 
develop  the  necessary  theory,  algorithms  and  computer  codes  for  simulating 
antennas  at  UHF  and  above  in  a complex  ship  environment.  A revised  mile- 
stone chart  for  the  three  year  proqram  is  shown  in  Table  I.  The  revised 
milestone  chart  shows  the  inclusion  of  the  more  useful  source  represen- 
tation litem  C of  the  Basic  Scattering  Code).  The  more  useful  source 
representation  is  being  given  priority  in  the  current  three-year  program 
over  the  near  field  analysis  of  coupled  structures  which  was  shown  in 
previous  milestone  charts. 

III.  BASIC  SCATTERING  CODE  DEVELOPMENT 

The  analytic  modeling  of  a ship  to  predict  the  radiation  patterns 
of  antennas  mounted  on  shipboard  is  being  accomplished  by  the  use  of  the 
Geometrical  Theory  of  Diffraction  IGTD)  [1,2,3].  This  is  a high  freguency 
technigue  that  allows  a complicated  structure  to  be  approximated  by  basic 
shapes  representing  canonical  problems  in  the  GTD.  These  shapes  include 
flat  and  curved  wedges  and  convex  curved  surfaces.  The  GTD  is  a ray  op- 
tical technique  and  it  theref--  ->  allows  one  to  gain  some  physical  insight 
into  the  various  scattering  and  diffraction  mechanisms  involved.  Conse- 
quently, one  is  able  to  quickly  seek  out  the  dominant  or  significant  scat- 
tering and  diffraction  mechanism  for  a given  geometrical  configuration. 

This,  in  turn,  leads  to  an  accurate  engineering  solution  to  practical 
antenna  problems.  For  example,  the  performance  of  an  antenna  at  UHF  and 
above  in  the  presence  of  a complex  ship  environment  can  be  modeled  using 
flat  plates,  box-like  structures  and  a finite  elliptic  cylinder  to  repre- 
sent the  various  component  parts  of  a ship  as  illustrated  in  Figure  1. 

A user  oriented  computer  code  is  being  developed  to  aid  the  antenna  and 
system  desiqner  to  accomplish  this  task. 

This  section  briefly  describes  the  basic  operation  of  the  third  gener- 
ation of  this  scattering  code  for  the  analysis  of  antennas  in  a shipboard 
environment.  This  version  of  the  code  allows  analysis  of  scattering  struc- 
tures that  can  be  modeled  by  concave  flat  plate  structures,  an  infinite 
qround  plane,  and  a finite  elliptic  cylinder  all  of  which  are  built  up 
from  the  basic  canonical  structures.  These  shapes  allow  one  to  model  a 
wide  variety  of  structures  in  the  UHF  range  and  above  where  the  scattering 
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TABLE  I 

REVISED  MILESTONE  CHART  (04/12/78) 


red  to  NOSC 


Sample  UHF  antenna  problem. 


structures  are  large  in  terms  of  a wavelength.  The  general  rule  is  that 
the  lower  frequency  limit  of  this  solution  is  dictated  by  the  spacings 
between  the  various  scattering  centers  and  their  overall  size.  In  practice 
this  means  that  the  smallest  dimensions  should  be  on  the  order  of  a wave- 
length. This  can  often  be  relaxed  to  approximately  a quarter-wavelenqth. 

The  GTD  approach  allows  other  types  of  solutions  to  be  interfaced 
with  it.  For  example,  the  modal  solution  can  be  used  for  a small  cylinder 
or  the  moment  method  can  be  used  to  model  wire  seqments,  small  structures 
or  antennas.  The  present  scattering  code  is  capable  of  interfacing  with 
the  wire  seqments  of  the  Numerical  Electromagnetic  Code  f NEC ) [4]. 

The  positive  time  convention  e'1(1)t  has  been  used  in  this  scattering 
code.  Also-, the  radiation  patterns  are  assumed  to  be  in  the  far  field. 

The  term  e”'1  /R  has  been  suppressed  unless  otherwise  specified  by  the 
user  of  the  code.  All  of  the  structures  in  this  code  are  assumed  to  be 
perfectly-conducting  and  the  surrounding  medium  is  free  space. 

As  mentioned  above,  the  GTD  approach  is  ideal  for  a general  high 
frequency  study  of  shipboard  antennas  in  that  only  the  most  basic  struc- 
tural features  of  an  otherwise  very  complicated  structure  need  to  be  mod- 
eled. This  is  because  ray  optical  techniques  are  used  to  determine  com- 
ponents of  the  field  incident  on  and  diffracted  by  various  structures. 
Components  of  the  diffracted  fields  are  found  using  the  GTD  solutions  in 
terms  of  the  individual  rays  which  are  summed  with  the  geometrical  optics 
terms  in  the  far  field.  The  rays  from  a given  scatterer  tend  to  interact 
with  other  structures  causing  various  higher-order  terms.  In  this  way 
one  can  trace  out  the  various  possible  combinations  of  rays  that  interact 
between  scatterers  and  determine  and  include  only  the  dominant  terms. 

Thus,  one  need  only  be  concerned  with  the  important  scattering  components 
and  neglect  all  other  hiqher-order  terms.  This  method  leads  to  accurate 
and  efficient  computer  codes  that  can  be  systematically  written  and  tested. 

Complex  problems  are  built  up  from  simpler  components  in  terms  of 
a modular  computer  code.  This  modular  approach  is  illustrated  in  the  block 
dianram  of  the  main  proqram  shown  in  Table  II.  The  scattering  code  is 
broken  up  into  many  subroutines  that  represent  different  scattered  field 
components,  ray  tracinq  methods,  and  shadowing  routines.  As  can  be  seen 
from  the  flow  chart,  the  code  is  structured  so  that  all  of  one  type  of 
scattered  field  is  computed  at  one  time  for  the  complete  pattern  cut  so 
that  the  amount  of  core  swappinq  is  minimized  thereby  reducing  overlayinq 
and  increasing  efficiency.  This  also  is  an  important  feature  that  allows 
the  code  to  be  used  on  small  computers  that  are  not  large  enough  to  accept 
the  entire  code  at  one  time.  The  code  can  be  broken  into  smaller  overlay 
seqments  which  will'  individually  fit  in  the  machine.  The  results  are, 
then,  superimposed  in  the  main  program  as  the  various  segments  are  exe- 
cuted. 


4 


The  code  is  divided  into  three  large  sections.  The  first  two  sections 
were  included  in  the  second  generation  code.  The  third  section  is  being 
included  in  the  third  generation  of  the  code.  The  first  section  contains 
the  major  scattered  fields  associated  with  the  individual  flat  plates  and 
the  interactions  between  the  different  plates.  These  include  the  direct 
field,  the  singly  reflected  fields,  doubly  reflected  fields,  the  singly 
diffracted  fields,  the  ref lected-diffracted  fields,  and  the  diffracted- 
reflected  fields.  The  diffracted  fields  include  the  normal  diffracted 
fields  as  well  as  slope  diffraction,  a newly  developed  heuristic  corner 
diffracted  field  and  slope-corner  diffracted  field.  The  double  diffracted 
fields  are  not  included  at  present,  but  a warning  is  provided  wherever 
this  field  component  might  be  important.  This  is  usually  only  a small 
angular  section  of  space.  This  field  may  be  included  later  whenever  time 
and  effort  permit.  The  second  section  contains  the  major  scattered  fields 
associated  with  the  elliptic  cylinder.  This  includes  the  direct  field, 
if  not  already  computed  in  the  plate  section,  the  reflected  field,  the 
transition  field,  the  deeD  shadow  fields,  the  reflected  field  from  the 
end  caps,  and  the  diffracted  field  from  the  end  cap  rims.  The  diffracted 
field  from  the  end  cap  rim  is  not  at  present  corrected  in  the  pseudo  caustic 
regions.  This  is  where  three  diffraction  points  on  the  rim  coalesce  into 
one.  This  is  only  important  in  small  angular  regions  in  space  and  is  not 
deemed  appropriate  to  be  included  at  the  present  time.  An  eguivalent  current 
method  could  be  used  for  this  small  region  but  it  is  rather  time  consumina 
to  use  for  the  benefits  derived  from  it  for  such  a general  code.  The  third 
section  contains  the  major  scattered  fields  associated  with  the  interactions 
between  the  plates  and  cylinder.  This  includes,  at  present,  the  fields 
reflected  from  the  plates  then  reflected  or  diffracted  from  the  cylinder, 
the  fields  reflected  from  the  cylinder  then  reflected  from  the  plates, 
the  fields  reflected  from  the  cylinder  then  diffracted  from  the  plates, 
and  the  fields  diffracted  from  the  plates  then  reflected  from  the  cylinder. 

These  terms  have  been  found  to  be  sufficient  for  engineering  purposes  when 
analyzing  wing-mounted  aircraft  antennas  as  well  as  many  other  structures  [5] . 

The  subroutines  for  each  of  the  scattered  field  components  are  all 
structured  in  the  same  basic  way.  First,  the  rav  path  is  traced  backward 
from  the  chosen  observation  direction  to  a particular  scatterer  and  subseguently 
to  the  source  using  either  the  laws  of  reflection  or  diffraction.  Each 
ray  path,  assuming  one  is  possible,  is  then  checked  to  see  if  it  is  shadowed 
by  any  structure  along  the  complete  ray  path.  If  it  is  shadowed  the  field 
is  not  computed  and  the  code  proceeds  to  the  next  scatterer  or  observation 
direction.  If  the  path  is  not  interrupted  the  scattered  field  is  computed 
using  the  appropriate  GTD  solutions.  The  fields  are  then  superimposed 
in  the  main  program.  This  shadowing  process  is  often  speeded  up  by  makina 
various  decisions  based  on  bounds  associated  with  '.he  geometry  of  the  structure. 
This  type  of  knowledqe  is  used  wherever  possible. 

The  shadowing  of  ravs  is  a very  important  part  of  the  GTD  scattering 
code.  It  is  obvious  that  this  approach  should  lead  to  various  discontinuities 
in  the  resulting  pattern.  However,  the  GTD  diffraction  coefficients  are 
designed  to  smooth  out  the  discontinuities  in  the  fields  such  that  a continuous 
field  is  obtained.  When  a scattered  field  is  not  included  in  the  result, 
therefore,  the  lack  of  its  presence  is  apparent.  This  can  be  used  to  advantge 
in  analyzing  complicated  problems.  Obviously  in  a complex  problem  not 
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all  the  possible  scattered  fields  can  be  included.  In  the  GTD  code  the 
importance  of  the  neglected  terms  are  determined  by  the  size  of  the  so- 
called  gliches  or  jumps  in  the  pattern  trace.  If  the  gliches  are  small 
no  additional  terms  are  needed  for  a good  engineering  solution.  If  the 
gliches  are  large  it  may  be  necessary  to  include  more  terms  in  the  solu- 
tion. In  any  case  the  user  has  a gauge  with  which  he  can  examine  the 
accuracy  of  the  results  and  is  not  falsely  led  into  believing  a result 
is  correct  when  in  fact  there  could  be  an  error.  In  order  to  show  the 
versatility  of  the  code  and  to  test  its  accuracy  a large  number  of  example 
problems  were  run  and  compared  with  measurements  whenever  possible.  These 
results  are  included  in  user's  manuals  [6,7]  that  illustrate  the  above 
ideas. 


The  method  used  to  input  data  into  the  computer  code  is  presently 
based  on  a command  word  system.  This  is  especially  convenient  when  more 
than  one  problem  is  to  be  analyzed  during  a computer  run.  The  code  stores 
the  previous  input  data  such  that  one  need  only  inout  that  data  which  needs 
to  be  changed  from  the  previous  execution.  Also,  there  is  a default  list 
of  data  so  for  any  given  problem  the  amount  of  data  that  needs  to  be  input 
has  been  shortened.  The  command  words  that  are  presently  imDlemented  in 
the  computer  code  are  listed  in  Table  III.  The  input  parameters  for  the 
third  generation  code  will  be  a modified  version  of  those  used  in  the  sec- 
ond generation  code  [7],  A more  detailed  description  of  the  modified  pa- 
rameters will  be  given  in  the  next  guarterly  report  and  the  user's  manual 
for  the  third  generation  code. 

The  input/output  format  of  the  basic  scattering  code  may  differ 
slightly  between  the  NOSC  and  OSU  versions.  NOSC  personnel  will  make  the 
modifications  in  the  input/output  for  their  needs  and  will  furnish  an 
appendix  to  the  user's  manual  detailing  their  modifications.  The  two 
versions  of  the  codes  will  be  the  same  other  than  the  differences  in  for- 
mat. 


In  order  to  confirm  the  validity  of  the  GTD  and  the  numerical  ac- 
curacy of  the  basic  scattering  code  over  a wide  range  of  geometries,  the 
code's  results  are  being  tested  against  measurements  made  at  OSU.  Whenever 
possible,  the  code's  results  are  also  being  checked  against  measurements 
made  at  other  facilities  or  against  results  found  in  the  literature.  For 
example,  the  basic  scattering  code  has  been  checked  against  measurements 
made  by  Bach  [8]  on  a variety  of  satellite  shapes.  The  code's  calculated 
results  for  a dipole  in  the  presence  of  an  eight  sided  box  as  illustrated 
in  Figure  2,  is  compared  against  measurements  in  Figure  3.  The  calculated 
results  for  various  locations  about  a finite  circular  cylinder  as  shown 
in  Figures  4a, b,  and  c,  are  compared  against  measurements  in  Figures  5a, b, 
and  c respectively.  In  all  cases  the  agreement  is  very  good,  confirming 
the  validity  of  the  basic  scattering  code's  results.  The  small  discrep- 
encies  that  exist  between  the  patterns  can  be  contributed  to  measurement 
inaccuracies  and  to  the  lack  of  double  diffraction  effects  in  both  the 
geometries  and,  in  the  case  of  the  cylinder,  to  the  pseudo  caustic  caused 
hv  the  coalescence  of  rays  on  the  endcap  rims. 
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TABLE  III 

COMMAND  WORDS  FOR  INPUT  DATA 

I 

CM:  Comment  card 

CE:  Last  or  only  comment  card 

TO:  Test  data  options 

UN:  Units  of  input 

FR:  Frequency 

PD:  Pattern  Data 

RG:  Far  Field  Range 

RT:  Translate  and/or  rotation  of 

geometry 

PG:  Plate  geometry  input 

GP:  Infinite  ground  plane 

CG:  Cylinder  geometry  input 

SG:  Source  geometry  input 

AM:  Source  input  from  NEC/AMP 

LP:  Line  printer  output 

PP:  Pen  plotter  output 

XQ:  Execute  program 

NX:  Next  problem 

EN:  End  program 


The  two  qeometries  considered  above  are  composed  of  either  plates 
alone  or  a finite  cylinder  bv  itself.  These  separate  shapes  can  be  an- 
alyzed usinq  the  second  qeneration  of  the  scatterinq  code.  The  third  qen- 
eration  of  the  code,  however,  has  the  pi ate-cvl inder  interaction  fields 
included  so  that  qeometries  containinq  both  plates  and  a finite  elliptic 
cvl inder  can  he  analyzed.  This  is  a much  more  difficult  problem,  due  to 
the  more  complicated  nature  of  equations  defininq  the  rav  paths  and  the 
fields.  In  order  to  test  the  accuracy  of  the  five  interaction  fields  that 
are  being  included,  a series  of  measurements  were  made  at  OSU  on  a simpli- 
fied qeometry  that  resembles  a situation  that  may  arise  on  shipboard. 

The  qeometry  is  shown  in  Figure  6.  This  is  a very  difficult  problem  be- 
cause of  the  multiple  bounces  of  the  fields  that  can  occur  between  the 
plate  and  the  cylinder.  This,  however,  illustrates  the  engineering  solu- 
tion that  can  be  obtained  using  the  scattering  code.  For  comparison  pur- 
poses, the  separate  pieces  are  first  studied  individually.  The  calculated 
and  measured  pattern  of  the  half-wavelength  dipole  is  compared  in  Figure 
7a.  The  results  for  the  square  plate  are  compared  aqainst  measurements 
in  Figure  7b.  The  results  for  the  finite  circular  cylinder  are  shown  in 
Figure  7c.  In  all  of  these  three  cases  the  comparisons  are  quite  good. 

The  radiation  pattern  of  the  dipole  in  the  presence  of  both  the  plate  and 
the  cylinder  is  compared  to  the  calculated  result  from  the  second  qener- 
ation scatterinq  code,  in  Figure  7d.  The  second  generation  code  does  not 
include  the  plate-cylinder  interaction  fields.  Although  the  basic  struc- 
ture of  the  pattern  is  there,  the  proper  levels  and  fine  details  are  mis- 
sing. The  measured  result  is  compared  to  the  third  generation  code  result 
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in  figure  7 e.  This  result  contains  five  pi ate-cvl inder  interaction  fields 
as  is  listed  above.  It  can  he  noted  that  the  levels  are  better  and  that 
the  fine  details  have  the  riqht  trends.  The  plate-cylinder  interations 
are  still  being  improved  and  more  tests  are  being  made  on  different  vari- 
eties of  shapes  to  see  if  the  code  is  functioning  properly.  However,  the 
results  in  Figure  7e  indicate  that  the  third  generation  code  is  givinq 
good  results  for  engineering  purposes. 

In  all  of  the  examples  above  the  scattering  code  uses  an  analytic 
representation  for  the  antenna.  The  antenna  presently  considered  in  the 
computer  code  is  an  electric  or  magnetic  radiator  with  a cosine  distri- 
bution, arbitrary  length  and  width,  arbitrary  magnitude  and  phase,  and 
arbitrary  orientation.  This  antenna  is  specified  as  follows  for  an  elec- 
tric source 


I(z* ) = I, 


■m 


jno  sinfl'  (‘fY  )cos(^  cosfi')  w sinfk^sinP'costb' ) -jkr' 

?tt  ^m  / .?  „ . oj  . n,  r1 

(jL-)  - cosV  kp-sinO'cosd)1 
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The  result  for  a magnetic  source  is  obtained  usina  duality.  Anv  arbitrary 
antpnna  can  be  simulated  bv  superposition  of  the  elements  by  makinq  the 
length  l small  < l 0.1X),  spacinq  the  elements  less  than  a quarter  wave- 
length, and  then  weighting  their  magnitudes  and  phases  to  simulate  the 
current  distribution  of  the  desired  antenna  [9],  Since  the  radiation 
pattern  is  relatively  insensitive  to  the  current  distribution  this  method 
works  very  well. 

The  NEC  code  is  also  interfaced  with  the  GTD  scattering  code.  This 
allows  complicated  antennas  to  be  analyzed  using  the  best  features  of  the 
moment  methods  to  find  the  antenna's  impedance  characteristics  and  current 
distribution.  Then  the  GTD  scattering  code  can  be  used  to  calculate  the 
radiation  pattern  of  the  antenna  in  the  presence  of  large  scattering  struc- 
tures. At  present,  the  scattering  code  is  written  to  use  only  the  wire 
segment  capabilities  of  the  NEC  code  and  not  the  patch  segments.  The  NEC 
currents  are  assumed  to  be  the  value  of  the  current  at  the  center  of  a 
segment  that  is  approximately  0.15X  or  less  in  length.  The  field  [ID] 
of  an  individual  segment  is  set  to 
jn0k«-  0- jkr ' 
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which  qives  consistent  results  with  the  NEC  code.  The  scattering  code 
is  normalized  in  the  far  field  to  exof-jkRl/R  unless  otherwise  specified. 

The  fields  are  output  in  volts/meter,  however,  it  should  be  noted  that 
unless  a range  R is  defined  the  units  are  actually  volts.  The  directive 
qain  is  also  calculated  by  using  the  power  radiated  as  supplied  by  the 
NEC  code. 

The  interfacing  of  the  NEC  and  GTD  codes  has  been  tested  by  using 
two  examples  obtained  from  Mr.  Jim  Loqan  at  NOSC.  Both  examples  are  composed 
of  antennas  that  were  modeled  by  wire  seqments  over  a perfectly-conducting 
infinite  qround  plane.  Using  the  assumption  that  the  current  distributions 
obtained  in  this  wav  will  not  be  perturbed  appreciably  by  slightly  chanqinq 
the  environment,  the  currents  from  the  NEC  code  can  be  used  to  simulate 
the  radiation  pattern  for  the  antennas  mounted  over  a finite  plate  which 
is  over  an  infinite  ground  plane  as  shown  in  Figure  8. 

To  test  the  compatibility  of  the  current  distribution  used  in  the 
GTD  code  with  the  current  weights  obtained  from  the  NEC  code,  the  GTD  code 
was  used  to  calculate  the  field  for  the  example  problem  of  four  quarter- 
wavelength  dipoles.  The  locations  of  the  four  dipoles  are  illustrated 
in  Figure  8 as  the  four  solid  lines  parallel  to  the  x-axis.  The  infinite 
ground  plane  is  coincident  with  the  square  plate  for  this  first  test, 
the  dipoles  were  divided  into  a total  of  20  segments.  The  fields  obtained 
from  this  calculation  resulted  in  a directive  gain  within  0.2  dB  of  the 
NEC  result  as  shown  in  Figures  9a  and  b for  the  plane  <t>=0°  and  <b=90°,  re- 
spectively. 

The  GTD  code  was  next  used  to  show  how  the  directive  gain  changes 
when  the  four  dipoles  are  placed  over  a square  plate  over  a qround  plane 
as  shown  in  Figure  8.  The  comparisons  between  the  infinite  ground  plane 
and  square  plate  cases  are  shown  in  Figures  10a  and  10b  for  the  two  pattern 
planes. 

The  second  example  is  composed  of  four  sets  of  crossed  dipoles  as 
shown  in  Fiqure  8 where  solid  and  dashed  lines  are  used  to  represent  the 
crossed  dipoles.  The  results  obtained  from  NOSC  for  this  example  were 
obtained  from  an  AMP  code  print  out  rather  than  an  NEC  code  print  out; 
however,  the  results  are  still  useful  for  comparison  purposes.  The  GTD 
and  AMP  results  for  48  segment  dipoles  over  an  infinite  qround  plane  are 
compared  in  Figures  11a  and  lib,  where  Figure  11a  shows  the  vertical  com- 
ponent and  Figure  lib  the  horizontal  component.  The  solid  curves  represent 
the  GTD  code  result  and  the  X's  show  the  AMP  result.  Again  the  two  results 
compare  within  0.2  dB.  Square  plate  results  are  also  shown  in  Figures 
11a  and  lib.  Directive  gains  for  the  infinite  ground  plane  and  square 
plate  cases  are  compared  in  Figures  12a  and  12b  for  the  <}>=0O  and  4>=45° 
planes,  respectively.  The  two  examples  confirm  that  the  GTD  and  NEC  code 
are  compatible. 
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Figure  R.  Geometry  for  the  Droblems  of  dipoles  over  a square 
plate  anH  infinite  qround  plane  showinq  the  side 
and  top  views. 
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Finure  9b.  Comparison  of  tho  directive  gain  of  four  diDoles  over  an 
infinite  around  plane  obtained  from  the  NEC  code  and  from 
the  GTD  code  using  the  NEC  currents  =90°  horizontal 
polarization) . 


Figure  10a.  Comparison  of  the  directive  gain  of  four  dipoles  over 
an  infinite  ground  with  four  dipoles  over  a square  plate  over 
a square  plate  over  an  infinite  ground  plane  (<t>  =0 
vertical  polarization) . 
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Figure  Ub.  Comparison  of  the  directive  qain  of  four  sets  of 
crossed  dinoles  over  an  infinite  around  and  over  a square 
elate  over  an  infinite  ground  4=0°,  horizontal  polarization) . 


Fiaure  i?a.  Comparison  of  the  total  directive  gain  of  four  sets  of 
crossed  diooles  over  an  infinite  around  and  over  a sauare 
nlate  over  an  infinite  around  f * =0° 'i . 
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F inure  1?h.  Comparison  of  the  total  directive  gain  of  four  sets  of 
crossed  dipoles  over  an  infinite  ground  and  over  a square 
plate  over  an  infinite  ground  (*=45°). 


The  main  disadvantage  of  using  the  NEC  currents  in  the  scattering 
code  is  the  time  required  to  compute  the  results.  At  present,  the  time 
for  one  pattern  computation  increases  bv  a linear  factor  of  the  number 
of  NEC  seqments  specified.  In  the  next  year,  ways  of  increasing  the  effic- 
iency of  this  orocess  will  be  studied.  A more  useful  source  representation 
capability  will  be  implemented.  The  inclusion  of  a width  to  the  source 
representation  is  a step  in  this  direction.  It  is  hoped  that  this  will 
be  compatible  with  the  patch  seqments  of  the  NEC  code.  Also,  a studv  will 
he  made  of  summing  up  the  incident  fields  from  each  of  the  source  dipole 
elements  so  that  the  diffraction  calculations  need  not  he  made  for  each 
individual  dipole.  The  effect  of  a qroup  of  dipole  sources  will,  thus, 
he  more  efficiently  computed. 

In  the  next  contract  year  the  third  generation  of  the  scattering 
code  will  be  completed  along  with  a user's  manual.  The  code  will  be  de- 
livered to  NOSC  by  October  31,  1978.  A code  manual,  which  has  already 
been  started,  will  be  written  giving  details  of  all  the  individual  sub- 
routines in  the  scattering  code.  Also,  as  mentioned  above,  a more  useful 
source  representation  will  be  studied  and  included  in  the  code  when  it 
is  completed. 

The  scattering  code  has  been  much  in  demand  over  the  past  year  and 
has  been  sent  to  many  agencies.  Because  of  its  versatility  it  can  be  used 
to  simulate  numerous  antenna  pattern  distortion  effects.  The  followina 
is  a partial  list  of  how  the  first  generation  code  has  been  used  and  by 
whom: 


Westinghouse  (baltimore,  Maryland)  is  using  the  code  to  analyze 
the  scattering  effects  of  a radar  (search  and  track)  mounted 
on  a tank. 

Electronic  Communications  Inc.  (St.  Petersburg,  Florida)  is 
using  the  code  to  determine  the  optimum  location  for  various 
antennas  mounted  on  and  around  a communication  van. 

Naval  Research  Laboratory  IWashinqton,  D.C.)  is  using  the  code 
to  locate  antennas  on  an  F - 1 4 aircraft. 

The  second  generation  code  has  been  or  is  in  the  process  of  beinq  sent 
to  the  fol lowinq: 

Martin  Marietta  Corp.  (Denver,  Colorado) 

Atlantic  Research  Corp.  (Rome,  New  York) 


IV.  REFLECTOR  ANTENNA  CODE  DEVELOPMENT 

The  task  of  the  reflector  antenna  code  development  is  to  develop 
user-oriented  computer  program  packages  by  which  the  fields  of  typical 
Navy  reflector  antennas  can  be  calculated.  These  codes  are  being  developed 
as  part  of  a larger  effort  to  develop  computer  models  for  simulating  an- 
tennas at  UHF  and  above  frequencies  in  a complex  ship  environment. 
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The  first,  task  tinder  the  reflector  antenna  code  development  was  to 
develop  the  computer  pronram  packaqe  for  calculatinq  far  field  patterns. 
Feed  hlorkuqe  and  sratterinq  effects  were  not  included  in  this  phase. 

These  will  he  included  later  t see  Table  I).  The  second  task  of  this  effort 
is  the  development  of  a second  qeneration  comnuter  code  hv  which  the  near 
field*  of  a typical  Naw  reflector  can  also  be  calculated.  The  second 
uenei  it  ion  reflector  code,  capable  of  both  near  field  and  far  field  compu- 
tations, will  be  completed  in  the  near  future  and  is  scheduled  for  delivery 
at  the  pnd  of  the  next  quarter  as  seen  in  Table  I. 

These  ctxiiDuter  proqram  packaqes  are  desiqned  to  have  the  flexibility 
to  compute  the  fields  of  most  reflector  antennas.  One  important  feature 
of  the  codes  is  their  capability  for  a general  reflector  rim  shape.  Another 
important  feature  is  their  capability  to  input  a practically  arbitrary 
volumetric  feed  pattern. 

Since  the  majority  of  Navy  reflector  antennas  have  parabolic  sur- 
faces, only  the  class  of  parabolic  surfaces  are  implemented  in  the  computer 
codes.  The  geometry  of  the  reflector  rim  is  treated  as  piece-wise  linear. 
The  codes  are  flexible  enough  to  include  general  reflector  rim  shapes  such 
as  ptical  and  rectangular  with  chopped  corners. 

.■  theoretical  approach  for  computing  the  far  field  pattern  of  the 
qeneral  reflector  is  haspd  on  a combination  of  the  Geometrical  Theory  of 
Diffraction  (GTD)  and  Aperture  Integration  (AI)  techniques.  AI  is  used 
to  compute  the  main  beam  and  npar  sidelobes;  GTD  is  used  to  compute  the 
wide-angle  sidelobes  and  the  backlobes.  To  implement  the  computer  algo- 
rithms based  on  these  theories,  efficient  ways  were  developed  to  handle 
calculations  involving  the  feed  pattern,  the  aperture  field  and  the  far 
fipld  pattern  computation. 

The  development  of  thp  first  generation  reflector  code  for  far  field 
patterns  was  completed  and  delivered  to  NOSC  in  February.  A user's  manual 
111]  has  hppn  published  for  the  first  aeneration  code.  Further  testing 
was' done  on  the  far  field  code  and  its  accuracy  has  been  verified  by  com- 
parisons with  independent  pattern  data.  The  results  of  these  comparisons 
are  qiven  later  in  this  report. 

A.  Feed  Patterns 

Sampled  data  from  each  measured  feed  pattern  cut  is  input  and  stored 
in  the  code.  Linear  interpolation  is  then  used  to  obtain  a piece-wise 
linear  representation  of  the  input  pattern  cut.  The  feed  patterns  in 
planes  other  than  those  corresponding  to  the  input  pattern  cuts  also  are 
calculated  by  linear  interpolation.  This  method  provides  a computationally 
efficient  way  of  calculating  the  aperture  fields  and  feed  spillover  radi- 
ation without  requiring  large  amounts  of  computer  storage  for  the  measured 
feed  pattern.  Only  relatively  few  data  points  need  to  be  stored  for  es- 
sentially complete  feed  pattern  information.  Furthermore,  the  piece-wise 
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linear  method  has  the  advantaqes  of  flexibility  and  simolicitv  for  qeneral 
feed  patterns.  No  cut-and-try  procedures  are  needed;  the  sample  feed  values 
can  be  obtained  directly  from  measured  feed  pattern  data. 

An  option  has  been  added  to  the  second  generation  reflector  code, 
currently  under  development,  for  convenience  in  specifying  the  feed  pat- 
terns for  offset-fed  reflector  antennas.  Part  of  the  loqic  for  the  offset 
feed  pattern  was  built  into  the  first  generation  code;  but  some  minor 
additions  were  necessary  to  run  the  offset  feed  option.  This  option  en- 
ables the  user  to  specify  the  patterns  of  the  feed  in  the  coordinate  system 
of  the  tilted  feed  horn  axis.  The  code  then  automatically  transforms  the 
feed  patterns  to  the  coordinate  system  of  the  reflector.  This  feature 
simplifies  the  input  of  data  for  a symmetrical  feed  horn  in  an  offset-fed 
reflector. 


Several  subroutines  have  been  developed  to  treat  a feed  with  a gen- 
eral pattern  and  general  polarization.  These  subroutines  decompose  the 
feed  radiation  into  components  to  provide  the  necessary  information  for 
the  following:  1)  aperture  fields,  ?)  GTD  edge  illuminations,  and  3)  direct 
feed  radiation.  Detailed  equations  for  the  feed  patterns  are  given  in 
the  Seventh  Ouarterly  Report  [IP]. 

The  far  field  reflector  code  uses  either  a piece-wise  linear  feed 
pattern,  or  the  analytic  functions  described  below.  For  sum  patterns  the 
analytic  function  is  given  bv 


where  the  constants  A,  4/  and  C can  be  controlled  for  each  input  pattern 
cut  <b  . The  pattern  valae  in  Equation  (1)  is  normalized  such  that 
F r (0 ) - 1 for  all  4>-plane  cuts.  The  constants  A,  C and  N control  the  shape 
of  the  pattern.  The  constant  permits  a given  pattern  shape  to  be  stretched 
or  compressed.  For  difference  feed  patterns  the  analytic  function  is  given 
by 


-A 

fA(4>)  = C e 
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B.  Reflector  Geometry  and  Aperture  Integration 

The  reflector  rim  geometry  is  specified  by  the  X-  and  Y-coordinates 
of  each  junction  point  on  the  rim.  Thus  the  projection  of  the  reflector 
rim  onto  the  aperture  plane  is  specified  as  piece-wise  linear.  The  rim 
points  P.  are  numbered  in  a counter-clockwise  sequence.  A typical  reflector 
rim  shape  with  eiqht  rim  points  is  shown  in  Fiqure  13a.  A rectanqular 
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grid  size  (D  and  D ) is  then  chosen  so  that  the  aperture  can  be  divided 
into  a principal  rectangular  grid  as  shown  in  Figure  13b. 

The  aperture  fields  are  calculated  and  stored  on  the  principal  grid 
for  use  in  the  aperture  integration.  The  principal  grid  values  are  used 
for  all  output  pattern  cuts.  The  aperture  fields  are  calculated  at  points 
off  the  principal  grid  by  using  linear  interpolation  from  the  principal 
grid.  This  is  more  efficient  than  calculating  the  aperture  fields  from 
the  feed  pattern  for  each  rotated  grid  that  is  used  for  off-principal  plane 
cuts. 

The  rotating  grid  method  described  in  Reference  [l3]  is  used  for 
the  aperture  integration.  This  method  provides  good  efficiency  because 
it  combines  several  features.  One  feature  is  the  use  of  overlapping  sub- 
apertures which  allows  a piece-wise  linear  representation  for  the  aperture 
distribution.  Thus  variations  in  the  aperture  fields  can  be  represented 
with  relatively  few  subapertures.  Furthermore,  the  subapertures  can  be 
electrically  large;  thus  minimizing  the  computer  storage  and  also  the  amount 
of  numerical  integration  required.  Another  feature  of  the  rotating  grid 
method  is  that  the  y-integrations  are  carried  out  for  each  column  of  the 
aperture  and  each  one-dimensional  integration  result  is  stored.  The  stored 
values  for  the  y-integration  are  then  used  for  each  pattern  angle  in  the 
plane  perpendicular  to  the  y-axis;  thus  the  efficiency  approaches  that 
of  a one-dimensional  integration.  Even  though  the  integration  grid  must 
be  rotated  to  obtain  the  pattern  in  other  planes,  the  required  grid  rotation 
is  computational ly  much  faster  than  the  numerous  two-dimensional  integrations 
that  would  otherwise  be  required. 

C.  GTD 

The  rotating  grid  method  is  used  to  calculate  the  far  field  pattern 
in  the  main  beam  and  near  sidelobe  region.  The  GTD  is  used  to  calculate 
the  wide-anqle  sidelobe  and  backlobe  regions  for  the  reflector  antenna 
computer  code.  Advantage  was  taken  of  the  existing  GTD  computer  subroutines 
used  in  the  basic  scattering  code.  The  GTD  analysis  of  the  reflector  is 
similar  to  that  of  diffraction  by  a flat  plate  [5],  except  that  the  curvature 
of  the  reflector  surface  must  be  taken  into  account.  It  was  found  that 
the  reflector  rim  must  be  subdivided  into  nearly  straight  segments.  A 
suitable  criterion  is  that  each  segment  of  the  reflector  rim  be  small 
enough  that  the  focus  lies  in  the  far  field  of  the  rim  segment. 

The  GTD  method  used  in  the  reflector  code  increments  around  the  rim 
and  determines  whether  a diffraction  occurs  for  each  linear  rim  segment. 

This  is  done  by  comparing  the  diffraction  angle  with  the  bounds  on  the 
permissible  range  of  angles.  If  no  diffraction  occurs  for  that  segment 
the  code  checks  the  next  rim  segment.  If  a diffraction  does  occur,  the 
diffraction  point  and  the  vector  for  the  incident  ray  from  the  feed  are 
calculated.  This  procedure  is  the  same  as  that  used  for  the  flat  plate 
scattering  code  except  that  the  geometry  information  associated  with  the 
parabolic  reflector  surface  is  changed. 
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Cach  diffracted  field  contribution  is  first  calculated  in  components 
perpendicular  and  parallel  to  the  specific  rim  segment.  The  perpendicular 
and  parallel  field  components,  EA  and  E„,  are  then  transformed  to  E and 
E components  in  the  reflector  coordinate  system,  so  that  the  diffracted 
field  contributions  can  be  summed  for  all  rim  seqments.  This  process  of 
calculating  and  summing  the  diffracted  fields  for  all  the  rim  seqments 
is  repeated  for  all  GTD  pattern  anqles. 

Most  of  the  far  field  pattern  calculations  are  performed  in  terms 
of  the  E and  E comoonents;  but  are  then  converted  to  principal  and  cross 
polarized  components  for  the  output  pattern.  A Huvqen's  source  represen- 
tation is  used  as  a reference  for  the  principal  and  cross  polarizationsflA] . 

D.  Summary  of  Computer  Code 

The  reflector  code  for  far  field  pattern  computations  requires  ap- 
proximately 1?F  Kbvtes  of  storaqe.  Typical  CPU  times  are  less  than  1 sec- 
ond per  pattern  anqle  on  the  ElectroScience  Laboratory  Datacraft  computer. 
For  example,  the  patterns  in  Fiqure  17  require  about  TOO  seconds  each  of 
CPI)  time  to  qive  a printout  for  36 1 pattern  values  (AA=0.5°).  The  cor- 
respondinq  CPU  time  for  Figure  23  is  about  210  seconds.  The  code  is  ex- 
pected to  run  3 to  5 times  faster  on  machines  comparable  to  a CDC-6F00. 

The  capabilities  of  the  code  may  be  summarized  as  follows: 

1.  A general  reflector  rim  shape  may  be  used  (piece-wise  linear). 

2.  The  required  input  data  for  the  feed  pattern  is  minimized  by 
piece-wise  linear  pattern  fittinq. 

3.  Storage  and  computation  time  of  aperture  data  for  AI  is  mini- 
mized by  using  a principal  rectangular  grid  and  interpolation 
of  the  aperture  field. 

a.  The  efficiency  of  the  far  field  pattern  computation  is  maximized 
by  the  use  of  GTD  for  wide  pattern  anoles  and  the  use  of  the 
rotatina  grid  method  for  near-axis  anqles  (main  beam  reoion). 

a.  The  feed  may  be  linearly  polarized  with  any  orientation  or 
circularly  polarized. 

The  GTD  and  Al  approaches  used  for  the  code  have  a basic  limitation 
on  the  minimum  size  reflector  that  can  be  modeled.  This  limitation  is 
nrohahly  on  the  order  of  1a  to  3x  for  the  reflector  diameter.  However, 
virtually  all  practical  reflector  antennas  exceed  3a  diameter.  There  is 
no  basic  limitation  on  the  maximum  size  of  the  reflector  for  the  basic 
analysis.  In  the  code,  the  reflector  surface  is  assumed  to  be  a perfect 
paraboloid.  Thus,  an  actual  reflector  antenna  must  have  sufficiently  qood 
tolerances,  especially  at  high  frequencies,  so  that  it  can  be  accurately 
modeled  by  the  code. 
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The  practical  limitations  on  this  version  of  the  code  can  he  sum- 
marized as  fol lows: 

1.  The  feed  must  he  located  at  the  focus  and  have  a constant  phase 
pattern. 

?.  Thp  reflector  surface  must  he  Daraholoidal . 

T.  Strut  scatterina  effects  are  not  included. 

4.  The  qrid  size  used  for  aperture  integration  must  be  chosen  suf- 
ficiently small  to  give  a good  representation  of  the  aperture 
field  distribution. 

f.  Arrav  variables  associated  with  the  rim  data,  the  principal 
qrid  and  the  feed  pattern  must  be  given  sufficient  dimensions 
for  the  required  input  data. 

Several  statements  are  included  in  the  code  to  print  out  a warning 
message  if  the  declared  dimension  is  exceeded  for  certain  arrav  variables. 
This  has  been  done  for  the  variables  mentioned  in  item  5 above. 

A block  diagram  of  the  main  program  is  shown  in  Table  IV.  A command 
word  system  is  used  to  input  data  into  the  code.  By  this  method  additional 
computer  runs  can  be  made  by  making  small  changes  in  the  input  data.  As 
seen  in  Table  IV  the  code  can  be  run  by  using  default  data  stored  in  the 
code  (NX:  Command).  This  is  useful  for  initial  testing  of  the  code.  Also, 
the  default  data  can  be  changed  as  desired  by  the  user  to  represent  a 
commonly  used  reflector  antenna.  Various  test  options  (TO:  Command)  are 
available  for  testinq  and  debugging  the  code.  The  input  patterns  of  the 
feed  antenna  are  specified  bv  the  desired  number  of  6-pattern  cuts.  The 
input  pattern  data  is  controlled  by  the  FD:  Command.  The  feed  pattern 
mav  he  specified  either  bv  sample  feed  data  points  or  as  analytic  func- 
tions. The  reflector  qeometrv  and  dimensions  are  specified  bv  the  DG: 
Command.  The  rectangular  qrid  size  to  be  used  for  aperture  integration 
is  also  specified  bv  the  DG:  Command.  A do  loop  for  operating  frequency 
was  included  so  that  the  code  could  be  run  for  different  frequencies  with- 
out chanqino  any  other  input  data  bv  using  the  FO:  Command.  This  is  useful 
if  the  feed  pattern  can  be  assumed  not  to  chanqe  with  frequency.  If  the 
feed  pattern  changes  with  frequency  the  feed  command  (FD:)  must  be  used 
before  the  code  is  executed  at  each  new  frequency.  The  6-planes  desired 
for  the  output  pattern  cuts  are  controlled  by  the  PH:  Command.  Either 
evenly  spaced  pattern  cuts  or  unevenly  spaced  cuts  (up  to  10)  can  be  spec- 
ified. The  LP:  and  PP:  Commands  provide  for  line  printer  and  pen  plotter 
outputs,  respectively.  Comments  can  be  printed  on  the  output  data  by  using 
the  CM:  Command.  The  UN:  Command  is  used  to  specify  the  appropriate  units 
for  the  input  data.  The  code  is  run  or  executed  by  using  the  XQ:  Command. 
The  command  word  system  is  described  in  detail  in  the  user's  manual  for 
the  code  [1 1]  . 
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TABLE  IV 

BLOCK  DIAGRAM  OF  THE  MAIN  PROGRAM  OF  THE  FAR  FIELD  CODE 


CM:  or  CErCommand 
[For  Comments  on  Output  ListingJ 


UN:Command 

Read  Units  f or  In put  Data 


XQ:Command 
! Execute  Program 


Write  Feed  Pattern 

If  1 LWFD) 

Calculate  Power  Radiated 
by  Feed  for  Antenna  Gain 

DO  Loop  on  Frequency 

I 

I 

I 


Calculate  Constants  Associated 
wi th  Basic  Reflect or  Geomet ry 


Set  up  Principal  Grid 

Calculate  and 
Fielc 

Store  Aperture 

Data 

Calculate  Geometry  Associated  with  GTD 
fGE0M  Subroutine! 

| DO  Loop  on  < 

i-plane  Cuts 

Calculate  Shadow  Boundaries 
( SBDY  Subroutine) 
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E.  Computed  Far  Field  Pattern  Results 


The  results  of  the  qeneral  reflector  code  for  several  examole  antenna 
qeometries  are  presented  in  this  section.  These  examples  illustrate  the 
versatility  of  the  code  for  computinq  far  field  patterns  of  reflector  an- 
tennas with  various  rim  shapes  and  feed  patterns. 

The  first  example  is  the  common  case  of  a circular  reflector  with 
an  on-axis  feed  as  shown  in  Fiqure  14.  A GTD  analysis  was  previously 
developed  for  this  case  as  reported  in  Reference  [15] . Calculated  results 
are  given  in  Reference  [15]  for  a 24"  diameter  reflector  with  an  F/D  ratio 
of  1/3  and  a frequency  of  11  GHz.  The  feed  is  linearly  polarized  in  the 
y-direction.  The  measured  feed  patterns  for  this  antenna  as  qiven  in 
Reference  [15]  are  shown  in  Figure  15.  The  piece-wise  linear  feed  pattern 
option  of  the  reflector  code  was  used  to  approximate  the  measured  H-  and 
E-plane  feed  patterns  as  shown  in  Figures  16a  and  b,  respectively.  The 
far  field  patterns  as  computed  by  the  general  reflector  code  are  shown 
in  Fiqures  17a  and  b for  the  H-  and  E-Dlanes,  respectively.  The  results 
from  the  general  reflector  code  were  found  to  be  in  good  agreement  with 
the  calculated  results  of  Reference  [15]  without  aperture  blockage  as  shown 
in  Fiqures  1R  and  19.  Aperture  blockage  and  feed  strut  scatterina  effects 
have  not  vet  been  included  in  the  qeneral  reflector  code.  The  comparison 
of  the  computed  results  with  the  measured  patterns  in  Fiqures  18  and  1Q 
is  one  indication  of  the  validity  of  the  qeneral  reflector  code.  Other 
comparisons  are  qiven  below  to  further  demonstrate  the  accuracy  of  the 
oeneral  reflector  code. 

The  next  two  examples  illustrate  the  flexibility  of  the  qeneral 
reflector  code  for  treating  various  rim  geometries.  One  example  is  a re- 
flector with  a 19X  x 1 1 A rectangular  aperture  and  a focal  distance  of 
1 OA . The  other  example  uses  the  19X  x 1 1 X rectanqular  aperture  but  with 
chopped  corners  as  shown  in  Fiqure  20.  The  feed  patterns  were  chosen  to 
be  the  same  for  both  reflectors  and  are  shown  in  Figures  21a  and  b.  The 
feed  is  vertically  polarized.  The  patterns  computed  by  the  general  re- 
flector code  are  shown  in  Figures  22a  to  c for  the  rectangular  aperture. 

The  corresponding  patterns  are  shown  in  Figures  23a  to  c for  the  same 
aperture  but  with  chopped  corners  as  shown  in  Figure  20.  Discontinuities 
occur  in  the  pattern  because  double  diffraction  effects  are  not  included. 
These  discontinuities  usually  occur  at  low  levels  in  the  pattern. 

It  was  originally  anticipated  that  offset-fed  reflectors  could  be 
handled  using  this  first  generation  code.  However  after  carefully  examininq 
this  code,  it  was  determined  that  several  improvements  should  be  made  before 
an  offset-fed  reflector  can  be  effectively  run.  Thus,  offset-fed  reflectors 
will  be  treated  in  the  second  generation  code  which  is  presently  beinq 
developed.  The  following  examples  for  offset  reflectors  were  run  with 
a preliminary  version  of  the  second  generation  code. 
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feed  Datterns  for  circular  reflector 
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Figure  17.  Far  field  patterns  of  circular  reflector  examole 
computed  by  general  reflector  code. 
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H-plane  pattern  of  a parabolic  reflector  with  a flanaed 
wavenuide  feed.  Computed  in  Reference  [15]. 
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Figure  22.  Far  field  patterns  of  rectangular  reflector  antenna. 
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Figure  22.  (Continued) 
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Figure  23.  Far  field  patterns  of  rectangular  reflector  antenna  with  chooped  corners. 
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Figure  23.  (Continued) 


The  first  example  for  an  offset  reflector  is  shown  in  Figure  24. 
Fortunately,  sufficient  data  were  published  in  Silver  [16]  to  completely 
specify  the  antenna  and  most  of  its  feed  pattern.  The  feed  pattern  data 
for  this  reflector  had  to  be  extracted  from  the  published  dB  contours 
[16]  of  the  aperture  field  as  shown  in  Fiqure  25.  The  feed  is  horizontall 
polarized.  The  feed  patterns  were  specified  by  using  the  option  for  an 
offset  feed  tilt  angle  (PSIT=20°)  and  the  feed  pattern  symmetry  option 
( I S YM=3 ) . This  feed  symmetry  option  specifies  the  feed  pattern  to  be 
symmetrical  about  the  y-z  plane.  Thus,  for  example,  the  feed  pattern  in 
the  $=-60°  cut  and  e>0  is  the  same  as  shown  for  6<0  in  the  *=+60°  cut  of 
Figure  26d.  Thus  the  input  feed  pattern  was  specified  in  the  $=-90  , 

-60  , -30°,  -15°,  +15°,  +30°,  +60°  and  +90°  $ plane  cuts  as  shown  in 
Figures  26a  through  e.  The  E-  and  H-plane  patterns  as  computed  for  9.2 
GHz  from  the  general  reflector  code  are  shown  in  Fiqures  27a  and  b.  The 
corresponding  measured  patterns  as  published  in  Silver  [16]  are  shown  in 
Fiqure  28.  the  measured  and  computed  patterns  are  in  good  aqreement  con- 
sidering that  feed  strut  scattering  effects  are  not  included  in  the  re- 
flector code.  The  full  180°  patterns  as  computed  by  the  general  reflector 
code  are  shown  in  Figures  29a  and  b.  The  feed  spillover  is  evident  from 
the  lobes  near  90°  and  110°  in  the  E-  and  H-planes,  respectively.  Since 
the  feed  pattern  information  was  obtained  from  the  aperture  field  contours 
of  Fiqure  25,  the  feed  spillover  was  estimated  by  extrapolating  the  feed 
patterns  with  smooth  curves  as  shown  in  Figure  26. 

The  next  example  is  the  offset  reflector  shown  in  Figure  30.  The 
projected  rim  shape  is  circular  as  shown,  althouqh  the  rim  is  actually 
elligtical.  The  feed  is  horizontally  polarized  with  an  offset  angle  of 
51.5  . A specialized  computer  code  was  previously  developed  by  Mentzer 
et  al  [17]  which  is  applicable  only  for  this  particular  antenna.  The  feed 
pattern  was  obtained  by  calculating  values  using  the  special  algorithms 
developed  in  Reference  [17];  and  then  these  values  were  input  using  the 
piecewise  linear  feed  option.  The  resulting  input  feed  pattern  cuts  are 
shown  in  Figures  31a  through  c.  Although  the  mid-C  band  frequency  was 
not  precisely  known,  several  frequencies  were  tried  until  one  was  obtained 
which  matched  the  beamwidths  of  the  patterns  calculated  in  Reference  [17]. 
The  patterns  computed  by  the  general  reflector  code  are  shown  in  Figures 
32a  to  c.  They  are  in  excellent  agreement  with  those  calculated  in  Ref- 
erence [17]  and  shown  in  Figures  33a  to  c.  An  antenna  gain  of  63  dB  was 
computed  by  the  general  reflector  code;  and  the  gain  of  the  feed  horn  was 
computed  to  be  about  37  dB.  Thus  the  general  reflector  code  predicts  the 
antenna  gain  as  26  dB  above  the  feed  horn  gain.  This  is  within  about  1 
dB  of  that  shown  in  Fiqure  33. 

The  depolarization  properties  of  offset  reflector  antennas  were 
studied  by  Chu  and  Turrin  [18].  They  calculated  the  maximum  cross  polari- 
zation produced  by  a linearly  polarized  feed  and  checked  their  result  with 
a measured  cross  polarized  pattern  of  the  offset  reflector  with  an  offset 
angle  of  A =45°  shown  in  Figure  34.  Kauffman  et  al  [lg]  also  checked  this 
case  with  calculated  patterns  from  their  computer  program  for  circular 


reflector  antennas.  Their  computer  program  is  based'  on  an  aperture  in- 
tegration analysis  involving  the  quantizing  and  sortinq  of  numerous  aper- 
ture field  data.  In  the  calculations  of  [18]  and  [19]  for  this  offset 
reflector  example,  a circularly  symmetric  feed  pattern  was  assumed.  The 
feed  pattern  was  taken  to  be  the  same  as  the  H-plane  pattern  of  an  open 
ended  circular  waveguide  as  shown  in  Figure  35.  The  patterns  computed 
by  the  general  reflector  code,  as  shown  in  Figures  36  through  38,  are  in 
excellent  agreement  with  those  computed  in  [1.9]  and  with  the  measured 
patterns  of  [18].  Note  that  excellent  aqreement  is  achieved  for  the  cross 
oolarized  beam  as  shown  in  Figure  37  for  the  4>=0o  plane.  There  is  no  cross 
polarized  beam  in  the  <j>=90°  plane  for  a symmetrical  feed  pattern. 


Figure  25.  Constant  intensity  contours  in  paraboloid 
aperture  (horn  feed  axis  titled  20°  relative 
to  paraboloid  axis). 

No  cross  polarization  occurs  when  the  offset  reflector  is  fed  bv 
a circularly  polarized  feed;  but  the  depolarization  property  manifests 
itself  bv  a beam  shift  for  circular  polarization.  The  left  and  right 
circularly  polarized  beams  shift  in  opposite  directions  [18].  The  cor- 
resnondinn  circularly  polarized  cases  as  computed  bv  the  qeneral  reflector 
code  are  shown  in  Fiqures  30  and  an.  The  calculated  beam  shift  in  each 
case  was  0.35°,  or  a total  shift  of  0.70°  between  left  and  riqht  circular 
polarizations.  This  compares  well  with  the  value  of  0.75°  as  determined 
from  measured  patterns  [18].  The  latter  are  reproduced  in  Figure  41. 
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la)  FRONT  VIEW  lb)  SIDE  VIEW 


Fiqure  30.  Offset  reflector  with  circular  projected  rim. 
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Fiaure  32.  (Continued) 


i Figure  33a.  4>  =0°  pattern  from  Reference  |7j. 

| 


83 


II 


(a)  FRONT  VIEW  (b)  SIDE  VIEW 


Figure  34.  Offset  reflector  geometry  for  example  of  Chu  and 
Turrin  [18].  D = 12  inches,  vi^ fi  =45°,  f/D=0.25 
and  freguency  = 18.5  GHz. 


Figure  35.  Feed  pattern  for  offset  reflector  example  of  Reference  [181. 
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Fiqure  36.  Far  field  pattern  for  example  of  Reference  [18], 
computed  by  general  reflector  code.  Principal 
polarization  is  <b=0°  plane. 
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Figure  37.  Far  field  pattern  for  example  of  Reference  [18], 
computed  by  general  reflector  code.  Cross  polarization 

in  d>=0°  plane. 
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Figure  38.  Far  field  pattern  for  example  of  Reference  [18],  computed 
bv  general  reflector  code.  Principal  Dolarization 
in  4=90°  plane. 
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Figure  39.  Far  field  pattern  for  example  of  Reference  [18], 
computed  by  general  reflector  code.  Left  circular 
polarization  in  4>=0°  plane. 
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Figure  40.  Far  field  pattern  for  example  of  Reference  p.R] , 
comDuted  by  general  reflector  cgde.  Right  circular 
polarization  in<t>=0°  plane. 
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Figure  41.  Measured  far  field  patterns  from  Chu  and  iurrin 
Circular  polarization  in  $=0°  plane. 
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V.  THEORETICAL  STUDIES 

Substantial  proqress  has  been  made  in  the  work  relatinq  to  the  ex- 
tension of  the  GTD  to  treat  the  diffraction  from  an  edqe  illuminated  bv 
a non  rav-optical  field.  In  the  conf iqurations  under  study,  the  non-rav 
optical  field  is  the  transition  reqion  field  of  a second  edqe.  This  topic 
has  been  described  earlier  on  pp.  67  and  68  of  [4608-6] ; it  has  also  been 
mentioned  on  p.  ?6  of  [4608-7]. 

It  should  he  pointed  out  that  Eauation  (41)  in  [4608-6]  applies  to 
a pair  of  staqqered  wedqes  with  different  wedqe  anqles,  as  well  as  to  the 
half  planes  shown  in  Fiqure  46  of  that  report.  In  the  present  period,  the 
analysis  has  been  extended  to  the  case  of  the  thick  edqe,  as  shown  in  Fiq- 
ure 4?. 


Fiqure  42.  Thick  edqe  illuminated  by  a source  at  0 or  O'. 


The  solutions  for  the  field  at  P when  the  source  0,  edqe  1,  edqe 
2 and  P are  all  in  line  can  be  summarized  as  follows. 

A.  The  electric  field  incident  from  0 is  perpendicular  to  the  top 
surface  (hard  boundary  condition). 

U(P)  = 1 ffc+d+r?)+nlk"1/?)  . 

B.  The  electric  field  incident  from  0 is  parallel  to  the  top  surface 
(soft  boundary  condition) 

U(P)  = l tan"!J  dTc+dVr?y  f(c+d+r?)  + Ofk-1)  . 

where  f(x)  is  defined  by  Equation  (42)  in  [4508-6]. 
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The  preceding  formulas  have  been  applied  to  an  example  where  the 
thick  edge  is  part  of  a rectangular  cylinder;  however  one  is  still  deal- 
inn  with  qrazinq  incidence  on  the  top  surface,  as  indicated  in  Fiqure  42. 
Calculations  of  the  field  at  P when  P is  in  the  far-zone  of  the  rectan- 
gular cylinder  are  comDarec!  with  those  obtained  from  the  moment  method. 

The  aqreement  is  found  to  be  very  close,  generally  within  a few  percent, 
even  when  the  surface  exposed  to  grazing  incidence  is  only  0.2A  wide. 

Also  the  solution  has  been  generalized  to  the  case  where  the  source 
(or  field  point)  is  not  in  line  with  edges  1 and  2,  as  in  the  case  of  the 
source  at  O' . 

This  research  is  essentially  complete  and  a technical  report  is  being 
prepared.  A paper  describing  this  work  was  presented  at  the  USNC/URSI 
Spring  Meeting  in  Washington,  D.C. 
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